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(57) Abstract: A method for the characterization of materials without using a physical model. In order to characterize a layer of 
material in a space A of values assumed by a function ?(?)(?: wavelength), (1) by reflectometry and/or ellipsometry on A, a measured 
spectrum ? is obtained, (2) m values ?i ? m of ? in A (m=l) are chosen, with B ={ ? such as min (? { ) =?= max (? d } for m>l, and B=A 
si m=I, (3) m complex index values n+jk to m ?i are chosen, (4) if m#l the index n (A) on B is calculated by interpolation, on the 
basis of (?j, n £ =n (?i)), l=i=m, and if m=l, n (?>=ni(? 1 ) on B, (5) M parameters M=2m-f 1 are chosen, in addition to an error function 
Er and by minimization of Er to M parameters, (a) with the aid of the law of interpolation of (?i,ni) on B, n (?), ? e B, (b) is deduced 
with the aid of n (?) and the thickness e of the layer, a theoretical spectrum ? (n (?) is calculated, e), (c) ? and ? are compared with 
the aid of Er and, ifi Er(?, ?) S e or minimal, it is possible to move on to (<u>e</u>), if not (d) M parameters are varied in order to 
tend towards the minimum of Er (?, ?) and it is possible to move on to (a), (e) if Er (?, ?)«u>e</u>, the index is equal to the latter 
thus obtained, if not m is increased and it is possible to move onto (2). 
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(57) Abrege : Procdde* de caracterisation optique de materiaux sans utilisation de modele physique. Pour caracteriser une couche 
d'un matenau dans un intervalle A de valeurs prises par une fonction a(X.) (A, : longueur d'onde), (1) par rSflectom&rie et/ou ellipso- 
m&rie sur A on obtient un spectre mesure" ^, (2) on choisit m valeurs cti ...o^de a dans A (m>l) , avec B ={ a tels que min (a {) <a< 
max (a {) } pour m>l, et B=A si m=l, (3) on choisit m valeurs d'indice complexe n+jk aux m <Xi, (4) si m#l on calcule par interpolation 
Tindice n (a) sur B, a partir des (<X|, nj=n (CO) . l^i^m, et si m=l, n (a)=ni(ai) sur B, (5) on choisit M parametres, M<2m+1, et une 
fonction d'erreur Er et par une minimisation de Er a M parametres, (a) a l'aide de la loi d'interpolation des (0ti,nj) sur B, on dSduit n 
(a) , a e B, (b) a 1' aide de n (a) et de Tepaisseur e de la couche, on calcule un spectre theorique ^ (n (a) , e) , (c) on compare ^ et * 
a 1' aide de Er et, si Er(^, $) S e ou minimal, on passe a (e), sinon (d) on fait varier les M parametres pour tendre vers le minimum 
de Er W) et on va a (a) , (e) si Er (*, *)<e, Tindice est pris 6gal a celui obtenu en dernier, sinon on augmente m et on va a (2). 
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PROCEDE DE CARACTERI S ATION OPTIQUE DE MATERIAUX SANS 
UTILISATION DE MODELE PHYSIQUE 

5 DESCRIPTION 

DOMAINS TECHNIQUE 

La pr^sente invention concerne un proc£de 
de caract^risation optique de materiaux. 

Ce proc£de permet de caracteriser des 
10 couches minces ou £paisses de ces mat^riaux, qui sont 
form^es sur des substrate. Les grandeurs physiques, que 
ce proc£d£ permet de determiner, sont: 

- l'ipaisseur d'une couche d f un mat£riau, 

- 1 1 indice de refraction de ce mater iau, et 
15 - le coefficient d 1 absorption de ce 

materiau. 

La caracterisation optique des materiaux 
est utile pour 1* analyse chimique de ces materiaux 
(notarament etude des bandes d f absorption, des 

20 proprietes de densif ication et des propriety 
d'oxydation) , dans les domaines de la micro- 
electronique, des capteurs, de la biologie, de la 
medecine) , ou bien pour 1' analyse des epaisseurs de 
d6p3ts de ces materiaux. 

25 Pour des exemples d ! application, on se 

reportera au document [1] qui, comme les autres 
documents cites par la suite, est mentionne a la fin de 
la prisente description. 

La caracterisation des proprietes optiques 

30 d'un materiau est egalement utile lorsque le materiau 
est par la suite structure ( pour y former par exemple 



WO 2004/061432 ^V/FR2003/050211 



2 

des gravures ou des rugosites) et que les proprieties de 
diffraction optique de la structure obtenue doivent 
etre calculees (voir le document [2] ) . 

Indiquons des a present que I 1 invention est 
5 particulierement utile lorsque la loi physique suivie 
par 1 1 indice de refraction complexe du materiau que 
l'on veut caracteriser n f est, a priori, pas connue. 



STAT DE LA TECHNIQUE ANTERIEURE 

10 On rappelle que les mesures optiques 

peuvent etre de diverses natures : 

II peut s'agir de mesures 

ref lectometriques • Dans ce cas, le coefficient de 
reflexion en intensity d'une structure est mesure sur 

15 un spectre (c f est~a-dire un intervalle) de longueurs 
d'ondes [A^, Am] . 

L 1 angle d 1 incidence de la lumiere 
d'eclairement peut §tre non nul. Le coefficient de 
reflexion peut etre mesure pour plusieurs angles 

20 d 1 incidence 0. Nous noterons R (0, X, p) le spectre 
ref lectometrique, od p est la polarisation du faisceau 
incident et X la longueur d'onde de ce dernier. 

Generalement , 1' angle 0 est nul et la 
polarisation p indeterminee . Dans le cas oft 0 n f est pas 

25 nul, il faut connaitre cette polarisation p. En 
general, cette derniere est de type (S) o€t (P) . 

II peut s'agir aussi de mesures 
ellipsometriques; Les grandeurs mesurees sont alors les 
parties reelle et imaginaire du rapport du coefficient 
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de reflexion en polarisation (P) au coefficient de 
reflexion en polarisation (S) . 

On note gen£ ra 1 ement p- I p I exp { j A) ce 
rapport coraplexe (avec j a =-l) . Les grandeurs 
5 gen§ralement traite"es sont ! p 1 , que l'on note tan(\jf) , 
et 009(A), ou des comb inai sons des deux. 

Par exemple, les variables issues d'un 
ellipsometre & modulation de phase sont les suivantes : 

I B =sin(2\|/) sin(A) et I c =cos (2\jf) . 
10 Un ellipsometre classique fournit, quant & 

lui, les variables suivantes : 

cMtan 2 \|f-l)/(tanV+l) ®t p^cosAd-oc 2 ) 372 . 

Par souci de g6neralit6, nous noterons Si et 
S 2 les grandeurs traities. 
15 Les spectres Si, ie [1,2], sont mesures sur 

une plage de longueurs d'ondes [A^, A«] . L ! angle 
d 1 incidence peut Stre quelconque. Plusieurs spectres 
peuvent itre mesuris £ diff brents angles d f incidence 
afin d'obtenir un spectre plus riche. Nous noterons 
20 s<9,A,)={Si(9,A,) ,S 2 (8A)} le spectre ellipsom^trique. 

De fa<?on complement a ire, des mesures 
goniom£triques (coefficient de reflexion en fonction de 
1' angle d f incidence) peuvent etre ajoute*es aux mesures 
servant & la caract^risation, afin de determiner 
25 l ! e*paisseur des di verses couches, pour une ou plusieurs 
longueurs d'ondes. Ces mesures ne sont pas suffisantes 
en elles-memes puisque l'on veut determiner 1 1 indice de 
refraction complexe sur une gatnrae spectrale allant de 
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Afin de siraplifier l 1 expose, nous noterons 
W un ensemble de spectres r£f lectomitrique ou/et 
ellipsometrique (s) (et eventuellement goniometriques 
pour quelques longueurs d'ondes) . 
5 Sans perdre aucune gen£ralite, nous 

n 1 exposerons , dans la present e description, le mode 
d'utilisation des proc6d6s de l'art ant^rieur et de la 
pr6sente invention que dans le cas d'une seule couche 
mince d 1 un mat€riau, f orraee sur un substrat connu. 
10 L^paisseur de cette couche est notee e et 

l f indice de refraction complexe du matiriau a la 
longueur d'onde X est not€ n* (A,) . 

On rappelle a ce propos que la partie 
r§elle (respectivement imaginaire) de cet indice de 
15 refraction complexe est notie n{X) (respectivement 
k(k)) et appel£e "indice de refraction" (respectivement 
"coefficient d' extinction" ) . 

En outre on note Er {*F (1) , ^? (2> ) une fonction 
d'erreur (par exemple l'6cart quadrat ique moyen) entre 
20 deux spectres *F (1) et ¥ (2> . 

Par exemple, on pourra prendre, lorsque 
l f on dispose de spectres ell ipsom^tr ique s sur plusieurs 
angles 9i, ie {l..n}, et d'un spectre r6f lectometrique : 

E r (Y (l) ,¥ l2> )o 

+ \R m (.Z)-R™(A)f\a (1) 



avec (X)={S X (1) (9iA) ,S2 (1) (6iA) ,R (l> (A.)} 
¥ ta) a)={S x <: " (6iA) ,S 2 «" (9 lf X) ,R ta > (X)} 
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et ie{l.„n} 

Des facteurs de ponderation peuvent £tre 
apportes a 1 1 integrale de f a<?on a ce que la f onction 
d'erreur puiBse tenir compte de variations sur la 
5 precision de mesure des spectres. 

La caracterisation optique de couches de 
niateriaux s f articule generalement autour de deux 
applications : 

La premiere application est le controle 
10 dimensionnel du depdt de couches minces que 1'on 
utilise en microelectronique . 

Generalement on connait bien le materiau 
depose, c'est-a-dire que l'on connait bien l'indice de 
refraction complexe de ce materiau aux longueurs 
15 d'ondes de la lumiere utilisee pour la caracterisation, 
Les lois suivies par 1 1 indice de refraction 
complexe sont soit tabulees, soit approchees par des 
lois physiques connues telles que, par exemple, le 
module de Cauchy, le modele de Sellmeier (voir le 
20 document [3] ) , les lois de Porouhi (voir le document 
[4]), et les lois lois d 1 oscillateurs harmoniques (voir 
le document [5] ) . Ces lois sont def inies par un notnbre 
fini de parametres. 

Par exemple, une loi de type Cauchy sans 
25 absorption, a deux parametres, est definie de la 
maniere suivante : 

Re[»*0l)] =n(X)=a 0 +^- 

Im[n*(A)]=k(X)«0 

Lorsque l f on est sur de la valeur des 
30 coefficients aids {0,1}) mais que I 1 on ne connait pas 
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l'epaisseur, un algorithrae de recherche est utilise 
afin de trouver l f 6paisseur qui minimise I'erreur entre 
la mesure Y et la r^ponse theorique Y compte tenu de 
1 1 indice modelis€. 
5 L'algorithme de recherche peut etre, par 

exemple, la raethode du Simplex, la recherche Tabou, la 
methode de Levendt- Mar quart ou la methode du recuit 
simuie (voir le chapitre 10 du document [6] ) . 

Lorsque 1 1 indice de refraction est 

10 approxiraatif , les coefficients ai sont integres dans la 
procedure d'ajustement de Y et *F . La recherche des 
coefficients a± constitue un procede de caracterisation 
de l f indice de refraction. 

Cependant, lorsque la loi suivie par cet 

15 indice de refraction est inconnue (il arrive que le 
materiau soit inconnu ou qu'il ne soit pas bien decrit 
par une loi physique connue) , ce procede reste 
approximatif et l'^paisseur risque d'etre fausse. 

La deuxi^me application est la 

20 caracterisation des materiaux. 

Le procede utilise reste le m§me, mis & 
part que le materiau n'est pas bien cormu. C'est 
precisement la fonction d 1 indice de refraction complexe 
la plus proche de la realite qui est visee. 

25 Le type de loi peut §tre choisi par 

analogie avec d'autres materiaux. Cependant , la loi 
suivie par 1' indice de refraction complexe peut §tre 
compliquee, ce qui est par exemple le cas d'une loi 
d'oscillateurs harmoniques : 

Ai Ai 

TTJL 

/=1 



30 k^)+/it(£)?=i+y 1 — — — — 

t^E + Ei + jTi E-Ei + jVi 
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Dans l f expression ci-dessus, j 2=s -l et 
l'indice de refraction et le coefficient d' extinction 
sont exprim£s non pas en fonction de X mais de E, avec 
E= 1240/A, (X en nm) . 
5 Dans ce cas, les coefficients des 

oscillateurs sont difficiles a trouver si l'on n l a pas 
leur ordre de grandeur. La recherche est dif f icilement 
automat isable, les algorithmes de recherche pouvant 
dormer des r£ponses erronees et le temps perdu peut 

10 §tre considerable. 

II existe une alternative a la recherche 
des coefficients : la methode point-a-point (PAP) . 
Cette methode PAP propose de ne pas choisir de loi 
physique et de rechercher l'indice de refraction 

15 complexe du materiau pour chaque longueur d'onde X±, ou 
ie [l...n], avec X 1 ^X l& et kn-Aw. 

Pour chaque A*/ un algorithme de recherche 
tente de trouver l'epaisseur, l'indice n(A,i) et le 
coefficient d' extinction k(A,i) qui minimisent l'erreur 

2 0 entre la raesure ¥ (A*) et la r^ponse theorique 

¥ (X if n(Xi) ,k(A,i) ,8) • 

Un tel proc^de pose un probldme parce que 
les divers points (Xi,E,n{Xi) ,k(Ad) ) ne sont pas 
forc^ment physiquement compatibles entre eux : par 
25 exemple, l'epaisseur trouvee peut varier en fonction de 
la longueur d'onde et la loi suivie par l'indice de 
refraction cotnplexe, plus simplement appeiee loi 
d'indice, peut presenter des discontinuites . 
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Ce proc^de est generalement valable 
seulement lorsque l'epaisseur est trls bien connue et 
que les mesures sont de trSs bonne quality. 

5 EXPOSE DE I.' INVENTION 

La pr£sente invention a pour but de 
reraedier aux inconvenients precedents . 

lie proc£d£ objet de 1' invention permet de 
caract€riser un mat£riau sans utiliser un module 
10 physique, C' est-&-dire sans utiliser une loi physique 
suivie par l'indice de refraction complexe du mat^riau 
e*tudi£. II est done tout particulidrement utile 
lorsqu'une telle loi n'est pas connue. 

Ce proc£d£ constitue une alternative aux 
15 proc£de"s de caract^risation connus, mentionnes plus 
haut. II peut Stre appele" "tn^thode des noeuds" car il 
utilise des "noeuds 11 c'est-£-dire des points de 
coordonn^es (X* n* ) , o\X n* est la valeur prise par 

l'indice de refraction complexe & la longueur d'onde X± 
20 et t prend un nombre limits de valeurs (enti§res) . 

De facon precise, la pre*sente invention a 
pour objet un # procid§ de caract#risation optique d'au 
moins une couche d'un materiau dans un intervalle A de 
valeurs prises par une fonction a d'une longueur d'onde 
25 optique X, lorsque X varie dans un intervalle de 
longueurs d'ondes, cette couche etant formed sur un 
substrat, ce proc^de" 6tant caract§ris6 en ce qu' il 
comprend les Stapes suivantes : 

1) on effectue un ensemble de mesures de 
30 r§f lectom6trie et/ou d» ell ipsom6 trie sur 1' intervalle 
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A, cet ensemble de mesures conduisant a un spectre 
mesur£, not£ et 1'on choisit les m^thodes de calcul 
associ^es a la nature des mesures et au type de couche 
a caract£riser, 

5 2) on choisit m valeurs initiales di -.0^ de 

la fonction a, appartenant a cet intervalle A, m 6tant 
un nombre entier au moins £gal a 1, et Ton d6finit un 
intervalle B comme £tant 1' ensemble des points a de 
1' intervalle allant du plus petit au plus grand des 

10 nombres 0d..-CX m/ lorsque m est sup£rieur a 1, et comme 
6tant 1' intervalle A lorsque m vaut 1, 

3) on choisit m valeurs initiales complexes 
d'un indice de refraction complexe n*=n+jk aux m points 
(Xi, i allant de 1 a m, 

15 4) lorsque m est different de 1, on choisit 

une loi d 1 interpolation qui permet de calculer l 1 indice 
de refraction n(a) du mat£riau sur 1' intervalle B, a 
partir des points (cc±, n±) , avec ni=n(<Xi), i allant de 1 
a m, et lorsque m vaut 1, n(a) est pris 6gal au nombre 

20 ni(<Xi) sur tout 1' intervalle B, 

5) on choisit M paramdtres variables, M 
itant infirieur ou §gal a 2m+l, 

6) on choisit une fonction d ! erreur 
Er(*F, x Sf) qui caract^rise l'6cart entre un spectre 

25 mesur§ Y et un spectre thiorique W , 

7) a l f aide d'une fonction de minimisation 
de Er^F, W) a M parametres, on effectue la s§rie 
d'£tapes suivante : 



WO 2004/061432 ^V/FR2003/050211 



a) a l'aide de la loi d' interpolation des 
(Oti.ni) sur 1' intervalle B, on d§duit n(a) , a 
appartenant a B, 

b) a l'aide de n(a> et de l'£paisseur e de 
5 la couche, et de m6thodes de calcul de spectres, on 

calcule un spectre thgorique W (n(a) ,e) , 

c) on compare *F et ¥ a l'aide de Er OF, W) 
et, si Er{*F, *F) est suf f isamment petite, c'est-a-dire 
inferieure a une valeur pr§d£finie e, ou est niinimale, 

10 on passe a l'6tape e) , sinon on passe a l'6tape d) , 

d) on fait varier les M param^tres 
variables de fa<?on & tehdre vers le minimum de Er(*F, 
*F) et l'on retoume a l'£tape a), 

e) si Er(Y, *F ) est inferieure a e, on 
15 obtient done un ensemble de M param&tres variables, 

pour lequel Er(*F, *P (n (0C,M) ,e) ) est minimum et l'indice 
de refraction est alors pris 6gal a celui qui a ete 
obtenu en dernier, et si Er(*F, *F) est sup^rieure ou 
§gale a e on passe a l'etape 8), 
20 8) on augment e le nombre m de valeurs 

initiales de la fonction a et l'on retoume a l'etape 
2) . 

II est done par exemple possible 
d'effectuer la caracterisation optique : 
25 - sur un intervalle de longueurs d'ondes X, 

en 1* occurrence 1' intervalle [Xmin, fanax] , 

- ou sur un intervalle de longueurs d'ondes 
inverses X/X, en 1' occurrence sur un intervalle 
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[(lA)min, (lA)max], ou (lA)min est Sgal £ 1/ (Jlmax) 
et (lA)max a 1/ (tanin) , 

ou sur un intervalle d' energies E (avec 
E = hv = he A oil h est la constante de Planck, c la 
5 vitesse de la lumidre dans le vide et V la frequence 
correspondant ^ X) , en 1' occurrence sur un intervalle 
[Emin, Emax] , ou Emin est §gal a he/ (tanax) et Emax a 
he/ (tanin) , 

- ou, plus g6n£ralement , sur un intervalle 
10 [ccmin, otmax] de valeurs prises par une fonction a de la 
variable X. 

II convient en outre de noter que 
1' invention est utilisable pour caract#riser un spectre 
ou une partie de spectre. 

15 Chaque loi d 1 interpolation peut etre 

choisie parmi les lois d 1 interpolation lin€aires, les 
lois d f interpolation cubiques, les lois d 1 interpolation 
polynomiales et les lois d 1 interpolation par exemple de 
type fonction « spline ». 

20 Selon un mode de mise en ceuvre prefer^ du 

procidi objet de 1' invention, le spectre est 
§chantillonn£ unif ornament en a (X) , e'est-a-dire que 
les valeurs initiales de la fonction a (voir l'^tape 2) 
mentionnSe plus haut) sont r^parties de manidre 

25 uniforme dans 1' intervalle A f la repartition des nceuds 
gtant ainsi homogdne. 

Corame on l'a vu, a{X) peut §tre choisi 
parmi X, 1/X et hc/X ou toute autre fonction de X, ou h 
est la constante de Planck et c est la vitesse de la 

30 lumi^re dans le vide. 
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De preference, a 1'etape 6) tnentionn^e plus 
haut, on mesure l'erreur sur un intervalle d' inter St C 
qui est inclus dans 1' intervalle B ou egal a cet 
intervalle B. 

5 lies M paramStres variables peuvent Stre les 

parties re*elles des indices de refraction aux points 
OLi, i allant de 1 a m, ou les parties imaginaires de 
ces indices de refraction, ou ces M param^tres 
variables peuvent Stre constitues par l'epaisseur du 

10 materiau dont on cherche 1' indice de refraction. 

La pre*sente invention a aussi pour objet un 
autre proce*de* de caracterisation optique d'au moins une 
couche d'un mater iau dans un intervalle de longueurs 
d»ondes [X rain, X max], cette couche etant formee sur 

15 un substrat, cet autre procede etant caracterise en ce 
que : 

- on effectue un ensemble de mesures de r£f lectometrie 
et/ou d'ellipsometrie, cet ensemble de mesures 
conduisant a un spectre mesure, note *F, 

20 - on choisit m longueurs d'ondes initiales Xi .JK* 
appartenant a cet intervalle, m etant un nombre entier 
au moins egal a 1, on associe, a chaque longueur 
d'onde, un indice de refraction, 

- on choisit une loi d ! interpolation au moins pour 
25 1 ! indice de refraction du materiau, pour les longueurs 

d'ondes comprises entre les longueurs d'ondes 
initiales X^Xa, 

- on choisit M param^tres initiaux, M etant au moins 
egal a m, a savoir un indice de refraction initial ni 

30 pour chaque longueur d'onde initiale A*, 1 ss i ss m, 
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les longueurs d'ondes initiales etant choisies de 
mani^re £. pouvoir determiner par interpolation au 
moins 1 ' indice de refraction pour toute longueur 
d'onde de l'intervalle [X min, X max], les couples (A*, 
5 ni) etant appeies noeuds, 

- on choisit des methodes de calcul de re"f lectometrie 
et d^llipsometrie, 

- on choisit aussi une fonction d'erreur Er, 
representative de I'ecart entre deux spectres *Fi et 

10 ¥2, les spectres *Fi et *F 2 etant calcules ou mesurls 
sur un nombre de points sup£rieur au nombre m de 
nceuds, 

- & I'aide des m longueurs d'ondes initiales, des M 
paramStres initiaux et de la loi d' interpolation , on 

15 met en ceuvre le processus d 1 optimisation suivant : 

- on determine un spectre theorique/ note *F , en 
fonction des methodes de calcul choisies, et de 
1' indice deduit par interpolation de sa valeur en A*, 
i allant de 1 a m, sur le spectre [A,min, Xmax] , 

20 — on determine l'erreur Er C¥, T ) , entre le spectre 
raesure et le spectre theorique, 

- on minimise cette erreur en faisant varier la 
position des valeurs des indices inconnus et/ou 
1'epaisseur de la couche et/ou les valeurs des indices 

25 de refraction aux longueurs d'ondes initiales, et l'on 
obtient un spectre, 

- on ajoute des longueurs d'ondes aux longueurs d'ondes 
initiales %x ~ ta, les longueurs d'ondes ajoutees 
const ituant de nouveaux noeuds, 
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— on repdte le procide en choisissant un nombre m' de 
longueurs d'ondes initiales, ra' etant superieur k m, 
et M' param£tres initiaux, M' etant superieur a M, 
jusqu'i. ce que la precision sur chaque spectre ainsi 
5 repr£sente au mieux soit egale 3. une precision 
predef inie. 

Dans ce cas, selon un premier mode de mise 
en oeuvre particulier, m est au moins egal & 2 ; selon 
un deuxierae mode de mise en oeuvre particulier, m est 

10 egal & 1 et l'on choisit des indices de refraction 
initiaux egaux. 

Dans ce cas egalement, selon un mode de 
realisation particulier, le materiau est non absorbant 
et le nombre M est egal k m, le coefficient 

15 d 1 extinction du materiau etant pris egal £. 0 ; selon un 
autre mode de realisation particulier, M est au moins 
egal & 2 m, on choisit en outre une loi d 1 interpolation 
pour le coefficient d' extinction du materiau, pour 
chaque longueur d'onde initiale Xi, 1 s i £ m, on 

20 choisit en outre un coefficient d 1 extinction initial 
ki, les longueurs d'ondes initiales etant en outre 
choisies de maniere & pouvoir determiner par 
interpolation le coefficient d 1 extinction pour toute 
longueur d'onde de l'intervalle [X min, X max] , et dans 

25 le processus d 1 optimisation, on minimise l'erreur en 
faisant varier en outre les valeurs des coefficients 
d' extinction aux longueurs d'ondes initiales et les 
longueurs d'ondes ajoutees sont en outre placees de 
facon £. representer au mieux le spectre du coefficient 

3 0 d' extinction du materiau. 
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Dans le cas de cet autre mode de 
realisation particulier, m peut §tre 6gal & 1 et l'on 
peut choisir des indices de refraction initiaux egaux 
et des coefficients d 1 extinction initiaux egaux. 
5 Encore dans le cas de cet autre proc6d6 

objet de 1' invent ion, la couche de materiau peut §tre 
mince, c'est-£-dire avoir une epaisseur inferieure a. la 
longueur de coherence de la lumi^re utilis^e pour les 
mesures, on peut choisir un paramdtre initial 

10 suppieraentaire, a savoir une Epaisseur initiale de 
couche, et dans le processus d» optimisation on peut 
minimiser l'erreur en faisant en outre varier la valeur 
de l r Epaisseur de couche ; dans une variante, la couche 
de materiau peut €tre epaisse, c'est-a-dire ne pas etre 

15 mince, et M peut §tre au plus egal a 2 m ; dans une 
autre variante, 1' epaisseur de la couche de materiau 
peut etre connue avec une precision suff isante et M est 
au plus egal a 2 m. 

La repartition des noeuds peut etre 

20 homog^ne. 

•-> 

BREVE DESCRIPTION DES DESSINS 

La presente invention sera mieux comprise a 
la lecture de la description d'exemples de realisation 
25 donnes ci-aprds, a titre purement indicatif et 
nullement limitatif/ en faisant reference aux dessins 
annexes, sur lesquels : 

- la figure 1 est une vue schematique de 
dispositifs permettant de caracteriser une couche 
30 conformement a 1' invention. 
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- la figure 2 montre les variations de 
l'indice de refraction en fonction de la longueur 
d'onde, pour un mater iau suivant une loi de Cauchy 
(courbe I) et pour un materiau caracterise conformement 

5 k 1' invention (courbe II) , 

- la figure 3 A (respectivement 3B> montre 
les variations de 1'indice de refraction 
(respectivement du coefficient d' extinction) en 
fonction de la longueur d'onde, pour un materiau 

10 suivant une loi a deux oscillateurs harmoniques (courbe 
I) et pour un materiau caracterise conformement 3. 
1 ' invention, et 

- la figure 4 illustre schimatiquement les 
paramdtres utilises dans une generalisation d'exemples 

15 de 1' invention . 

EXPOSE DfiTAILIiS DE MODES DE REALISATION PARTICULARS 

L 1 invention propose une alternative aux 
procedes classiques, mentionnes plus haut. Elle permet 
20 d'allier la coherence d'un module de couche 
(correspondant £ une loi d f indice continue et £ des 
epaisseurs physiques constantes) , Zk la generalite 
concernant la loi d'indice a trouver (comme dans la 
methode PAP) . En outre, la resolution n'est limitee que 
25 par la resolution du spectre mesure. 

Dans le procede objet de l f invention, le 
spectre d'indice n* (X) est caracterise par : 

- un nombre reduit de "nceuds" , qui sont des 
points de coordonnees (Xi,ni,ki) ou (Xi, n*i) , avec ni = n 

30 (%0 , ki a k (^i) et n*i=ni+jki, ou j 2 «-l, et 
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- une loi d 1 interpolation entre les nceuds, 
qui peut §tre, par exemple, lin€aire, cubique, de type 
« spline » ou polyndmiale (de degr§ quelconque) . 

Cette loi d 1 interpolation permet de 
5 calculer, h partir des nceuds, les indices de refraction 
et les coefficients d 1 extinction pour les longueurs 
d^ndes situ^es entre les nceuds. 

Par exemple, lorsque 1 1 indice de refraction 
est caract£ris6 par un ensemble de valeurs pour des 
10 longueurs d'ondes %i ~ X^, on peut utiliser une 
interpolation lin£aire entre deux longueurs d f ondes Xi 
et X±+ x pour calculer 1' indice n & la longueur d'onde X 
(voir le document [6] chapitre 3) : 

15 avec X±<X<Xt+i 

On peut faire de mSme pour le coefficient 
d' extinction. 

Lorsque le nombre de nceuds le permet, des 
formules d 1 interpolation plus complexes, faisant 
20 intervenir les nceuds voisins, peuvent §tre utilis£es 
(voir le document [6] chapitre 3) . 

Un module de couche est done caract§ris£ 
par une epaisseur e et une famille de nceuds. 

On donne ci-apres un exemple du proc§de 
25 objet de 1' invention. Dans cet exemple, les raesures *F 
sont constitutes d'une mesure r6f lectom^trique R(X) et 
d'une mesure ellipsom€trique S 1#2 O/W oil 0 est 1" angle 
d» incidence du faisceau lumineux que I 1 on envoie sur la 
couche & etudier lors de la mesure d 1 ellipsom6trie. 
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Cette couche est une couche mince de sorte 
que l'gpaisseur de cette couche est aussi une variable 
du probleme. En outre, on suppose qu'une seule couche 
est inconnue, cette couche 6tant form^e sur un substrat 
5 connu . 

Expliquons d'abord briSvement cet exemple, 
qui utilise un algorithme (algorithme de la < m^thode 
des noeuds » conformlment a 1' invention) > 

A partir d 1 informations suppos€es sur 
10 1 1 £paisseur e de la couche 6tudi£e et sur 1 1 indice de 
refraction n(X) et le coefficient d f extinction k(A,) du 
mat^riau de cette couche, on construit des noeuds de 
depart (en faible nombre) et une £paisseur de depart e. 

On dispose ainsi de m noeuds et , par 
15 interpolation, on peut connaitre n(X) et k(k) en dehors 
des valeurs des longueurs d'ondes associ£es aux noeuds. 

A partir de l'€paisseur de depart 8 et de 
ces valeurs de depart n(X) et k(A,) , on determine le 
spectre th£orique W en utilisant des calculs 
20 ellipsomitriques et r£f lectrom^triques . 

Par ailleurs, au raoyen de dispositifs 
d'ellipsometrie et de r§f lectrometrie et d'un 
spectromdtre, on obtient Si, 2 (6 ,X) et R(X) et I'on en 
diduit les mesures noties ¥ (pour des conditions de 
25 mesure 8 et 31) . 

On compare ensuite ¥ et f en utilisant une 
fonction d'erreur Er et l f on optimise la valeur de 
1 1 indice de refraction et la valeur du coefficient 
d 1 extinction aux diff^rents noeuds, ainsi que la valeur 
3 0 de I'gpaisseur, en cherchant a minimiser Er (*F 5 \F) . 
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Lorsque ces valeurs sont optimisees et si 
la precision sur le spectre n(k), le spectre k(k) et 
l»epaisseur £ n'est pas suffisante / on ajoute de 
nouveaux noeuds, on fait varier l'epaisseur 8, et 1'on 
5 recommence la determination de ¥ , la comparaison de 

et *F et 1 •optimisation que 1'on a mentionnees plus 
haut , etc . 

On arrite la boucle ainsi d£finie lorsque 
la precision sur chacun des spectres n(X) et k(A,) et 
10 sur l'epaisseur 8 est jug6e suffisante (ajustement 
satisfaisant de ¥ et . 

Les spectres n(X) , k{X) et l'epaisseur e 
sont ainsi caract6ris£s. 

Dans un exemple donne a titre purement 
15 indicatif et nullement limitatif, on trouve une 
epaisseur £ egale a 212, 3nm pour la couche. 

La figure 1 montre de fa<?on schematique la 
couche etudiee 2, formee sur un substrat 4. On voit le 
dispositif d f ellipsom£trie 5,6, le dispositif de 
20 r£f lectrom€trie 8 et le spectrometre 10. On voit en 
outre des moyens £lectroniques de traitement 12, 
comprenant tin ordinateur et permettant de caracteriser 
n(A,), k(A-) et £ en fonction des informations fournies 
par le spectromStre 10 et conformement au proc§d£ de 
25 1 1 invent ion . 

Ces moyens 12 sont munis de moyens 
d'affichage 14 qui permettent, en particulier, 
d'afficher la courbe des variations de n en fonction de 
X et la courbe des variations de k en fonction de X. 
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Revenons de facon plus detaillee sur 
1'exemple donne. 

Phase 1 

lie proced£ de cet exemple comporte d ! abord 
5 une £tape d ! initialisation. 

L 1 algorithme debute avec un nombre reduit 
de noeuds, plus precisement au moins un noeud. On peut 
done debuter par un seul noeud, en iraposant un indice de 
refraction et un coefficient d 1 extinction qui restent 
10 constants lorsque la longueur d f onde varie. 

On choisit des positions de noeuds de facon 
i. pouvoir, k partir de cette famille de noeuds, diduire 
tout le spectre par interpolation. Le module de couche 
est done k 3 pararaetres ou plus, puisque l'epaisseur 
15 est aussi une variable i. determiner. La table d 1 indice 
sur tout le spectre est done d€duite des noeuds par 
interpolation. 

C'est le cas lorsque l'epaisseur des 
couches est de 1 1 ordre de grandeur de la longueur 
20 d'onde (couches minces) . Mais lorsque l'epaisseur de la 
couche est sup^rieure eI la longueur de coherence de la 
source de lumiere, l'epaisseur n 1 intervient quasiment 
plus dans le calcul de la r6ponse de la couche, et 
n'est done plus une variable du probleme. 
25 Par exemple, dans le cas de disques 

optiques (CDROM) , sur lesquels une couche tres 
importante est d6pos£e (son epaisseur est de 1' ordre de 
1 millimetre), le coefficient de reflexion d'une telle 
couche n'est plus fonction de I'epaisseur de la couche 
30 mais seulement de I 1 indice de refraction de celle-ci, 
la longueur de coherence du faisceau lumineux incident 
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6tant inf£rieure a l'6paisseur de la couche. Dans ce 
cas precis, la longueur de coherence du faisceau 
incident est d£termin£e par la rugosite des couches. 

On choisit/ par exemple, de placer les deux 
5 premiers noeuds aux extr€mit£s X min et X max du 
spectre. Les valeurs de l f indice complexe en ces 
extr6mit£s sont choisies en fonction du type du 
materiau £tudi£. Par exemple, sur un spectre 
ellipsom^trique entre 300nm et 800ntn d'une couche mince 
10 de r£sine photosensible, on prend n(300nm) =n(800nm) =1, 5 
et k(300nm) =k(800nm) «0. 

Lorsque le spectre n'est caract#ris£ que 
par deux noeuds, l'indice entre les extremes est 
determine par interpolation lin^aire. Dans le cas 
15 consid6r6, on a done n(X)=l,5 et k(^)=0 pour X 
appartenant a [300nm, 800nm] . 

A partir de trois noeuds, on choisit plut6t 
une interpolation cubique afin d'obtenir des formes de 
loi d'indice plus douces que celles obtenues par 
20 interpolation lin6aire. 

L'epaisseur de depart est, quant a elle, 
choisie aussi proche que possible de l'6paisseur 
r6elle. 

Phase 2 

25 On proedde ensuite a une determination 

optimale des valeurs de l»indice de refraction et du 
coefficient d 1 extinction sur les noeuds et de la valeur 
de 1 1 ^paisseur . 

Pour ce faire, les spectres *P (X) Bont 
30 calculus a I'aide du module de couche utilise, 
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resultant du choix des noeuds, de la loi d 1 interpolation 
et de l f epaisseur de la couche. 

Le moddle physique utilise pour le calcul 

de *F est bien stir fonction de la mtthode de mesure 
5 utilisee, c'est-S.-dire notamment de l 1 angle d 1 incidence 
de la luraidre, du spectre utilise et du module de 
couches minces ou de couches tpaisses le cas 6ch6ant 
(voir par exemple les modules de couches empil^es dans 
le document [3] ) . 

10 Les spectres ¥ 6tant constitute par un 

ensemble de mesures de natures diverses (par exemple 
des mesures ellipsomttriques et r§f lectomttriques) , on 
utilise pour les mesures de ref lectomitrie 
(respectivement d 1 ellipsom£trie) une mtthode de calcul 

15 de r£f lectometrie (respectivement d 1 ellipsomttrie) • 

Les mesures r£f lectometriques et 

ellipsom^triques sont combinees par 1 1 intermtdiaire 
d'une fonction d f erreur ErtT,^) qui est par exemple du 
genre de celle qui est d^finie par liquation (1) . 

20 Gri.ce k une fonction de recherche, on 

minimise I'ecart entre ¥ (X) et ¥(A.) , en faisant varier 
la valeur de 1 1 indice de refraction et la valeur du 
coefficient d 1 extinction k la position de chacun des 
noeuds ainsi que l'tpaisseur de couche (si cette 

25 £paisseur est un facteur important dans le calcul de 
Lorsque I'ecart est minimum c'est-^-dire lorsque 
Er(*F,*F) est minimum/ cela signifie que les mesures 
r£flectom£triques et ellipsomttriques coincident au 
mieux (pour un nombre de noeuds donne) . 
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A ce stade, on obtient, pour un nombre 
connu de nceuds et une position spectrale connue pour 
chacun de ces des nceuds, le module de couche (indice de 
refraction, coefficient d' extinction et epaisseur) qui 
correspond au mieux k la couche reelle. 

La validity du module trouv6 est d'autant 
plus assured que le nombre de mesures est grand. Pour 
avoir un grand nombre de mesures, on peut par exemple 
utiliser plusieurs angles d 1 incidence de la lumidre 6i # 
1 as i & l t pour les spectres ellipsometriques, faire 
une mesure ref lectometrique et faire des mesures 
goniometriques complementaires . 

Phase 3 

Ensuite, on augmente le nombre de nceuds. 

On ajoute un nombre fini de nceuds. Dans un 
premier mode de realisation, les nceuds rajoutes sont 
positionnis de fagon h representer au mieux les 
spectres n{X) et k(A,) . A titre d 1 exemple, on place ces 
nceuds supplementaires aux endroits od l'^cart entre *P 

et Y est maximum ou aux endroits en lesquels les nceuds 
sont les plus espac^s. Et I'on retourne & la phase 2 
tant que la precision sur chacun des spectres n(A,) et 
k(k) et sur 1' Epaisseur e n'est pas suffisante, 
c f est-&-dire n'est pas egale a une precision 
predef inie. 

Dans un deuxi&me mode de realisation, on 
intercalle de nouveaux nceuds entre deux nceuds de 
1' ensemble de nceuds pr^alablement choisis, ces nouveaux 
noeuds se repartissant de fa<?on homogene sur le 
spectre . 

II convient toutefois de noter ce qui suit. 
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Lorsque l'on augmente le norabre de noeuds, 
il n'est pas necessaire de conserver la position des 
anciens nceuds. Par exemple, adraettons que l'on 
6chantillonne unif ormiment un spectre de 400nm a 800nm 
5 avec 3 noeuds, Ces noeuds se situent done respectivement 
a 400nm / 600nm et 800nm. Lorsque l'on passe a 6 noeuds, 
les trois noeuds additionnels ne peuvent etre places de 
facon a ce que le spectre soit unif ormement 
ecliantillonn£ si l'on conserve la position dee anciens 

10 noeuds. La position des 6 noeuds peut £tre definie, si 
l'on souhaite un echantillonnage uniforme par les 
valeurs 400, 480, 560, 640, 720, 800nm. L'ancien noeud 
central a 600nm disparait done. Pour calculer la valeur 
de l'indice a ces positions a partir des anciens noeuds, 

15 on procede par interpolation. 

Dans ce qui suit, nous donnons deux 
exemples courants d» application de l 1 invention. Ces 
exemples mettent en jeu deux types de materiaux qui 
suivent des lois differentes. 

20 A partir de raesures ellipsometriques et 

reflectom§triques, nous proposons de retrouver les lois 
physiques suivies par ces materiaux. 

Nous proc€dons de la facon suivante : 

Un mater iau fictif est cre6, ce materiau 

25 suivant une loi theorique connue (une loi de Cauchy ou 
une loi d'oscillateurs harmoniques) , avec des 
parametres que nous fixons arbitrairement . Les 
variations de 1 1 indice de refraction coraplexe en 
fonction de la longueur d'onde sont ainsi parf aitement 

30 connues. 
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Nous imposons en outre au materiau une 
6paisseur de 200,00nm, sur un substrat de silicium, ce 
dernier etant egalement tr^s bien connu. 

Des raesures fictives (mesures 

5 ellipsometriques, mesures ref lectometriques) sont 
calculees, puis bruitees de facon & introduire un 
defaut d*appareillage. 

Tout se passe comme si nous avions des 
mesures reelles, ef features sur le materiau. Mais, 
10 contrairement k la realite, nous connaissons 
parfaitement l'indice de refraction complexe puisque 
nous l'avons fixe, de meme que nous avons fixe 
l'epaisseur de la couche de materiau. 

Nous testons ici la methode "en aveugle", 
15 c'est-H-dire que nous partons d ! une epaisseur . fausse 
(220nm) et d» indices de refraction complexes faux, 
puisqu*ils sont senses £tre inconnus. 

Nous appliquons le procede objet de 
1' invention puis nous coraparons l ! indice de refraction 
20 complexe trouve k l'indice de refraction complexe 
th^orique. Nous retrouvons bien les mSmes lois, tres 
precis^ment, ainsi que la m£me epaisseur de couche. 

Comme premier exemple, prenons un materiau 
dont I'indice de refraction complexe suit une loi de 
25 Cauchy telle que : 

/<lx ^ rt ,300 2 _300 4 
k(A,)=0 

Cette loi d'indice est typique des resines 
photosensibles (dans la plage spectrale allant de 300nm 
30 k 800nm) . 
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Afirx de retrouver £ l ! aide de la methode 
des nceuds (c 1 est-^-dire du procede objet de 
l 1 invention) la loi d'indice mentionnee ci-dessus, nous 
effectuons deux mesures, a savoir une me sure 
5 ellipsom£trique a un angle de 70° et une me sure 
r£f lectometrique . 

Les conditions de traitement sont les 

suivantes : 

- le spectre traits est compris entre 300nm 

10 et 800nm, 

- au depart, les nceuds sont aux positions 
(300nm, 1,6) et (800nm, 1,6) , c'est-^-dire que l f indice 
est considere comme variant lin^airement entre 300nm et 
800nm, et que sa valeur est constante (egale k 1/6), 

15 - le nombre de nceuds est augmente 

it^rativement selon la sequence 2— »4— >6, 

- lors de la procedure d' augmentation des 
nceuds, la position en longueur d'onde de chaque nceud 
est calculie de fa<?on que 1'^chantillonage en 1/X soit 

20 uniforme (X : longueur d'onde) , la densite de points 
£tant done d' augmented vers les faibles longueurs 
d'ondes, 

- la loi d' interpolation est une loi 
cubique, lorsque le nombre de nceuds est superieur 3. 2, 

25 sinon elle est liniaire, et 

- l'algorithme de minimisation utilise est 
un algorithms de type Simplex. 

La figure 2 permet de comparer 1 ' indice de 
refraction correspondant au mate"riau fictif qui suit 
30 parfaitement une loi de Cauchy (courbe I) Bl 1' indice de 
refraction que nous trouvons par la methode des nceuds 
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(courbe II) , k l'aide de 6 noeuds (representee par" des 
cercles sur la figure 2) . Nous avons utilise un *F 
compost d'une mesure ellipsometrique {Si {X) , S 2 (X) } a 70° 
et d'une mesure ref lectometrique R(X) . 
5 Nous aboutissons & une epaisseur de 

199,838411m. 

Comme deuxieme exemple, prenons un materiau 
dont l'indice de refraction complexe suit une loi £. 
deux oscillateurs harmoniques, telle que : 

J i-tE + E t +jG t E-E,+JG t 

avec 

j 2 = -1 et E = 1240/X (X en nm) 
Ai=0 , 25X1 , 5 2 A 2 =0 , 25X1 , 5 2 
E!=1240/400 E 2 «1240/300 
15 Gi=0 f 3 G2=0,3 

Dans ce deuxieme exemple, la methode des 
noeuds est appliquee a un ensemble de mesures 
ellipsometriques effectuees entre 250nm et 800nm, k 
75°, 70° , 60° et 45°, avec en plus une mesure 
20 ref lectometrique. L« epaisseur reelle du materiau etant 
de 200nm, on trouve une epaisseur de 200,25 nm avec la 
methode des noeuds. L'ajustement sur la loi d'indice 
consideree danB ce deuxidrae exemple est tres bon, comme 
le montrent les figures 3A et 3B. 
25 Ces figures 3A et 3B illustrent 

respect ivement les reconstructions des courbes n(X) et 
k(X) du materiau i. l'aide de la methode des noeuds. La 
reconstruction est effectuee k I'aide de quatre 
spectres ellipsometriques et d'un spectre 

30 ref lectometrique . Les pics d ! absorption reels sont tr£s 
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bien repr€sent£s par la courbe obtenue par 
interpolation cubique entre les noeuds (represents par 
des cercles sur les figures 3A et 3B) . 

Sur la figure 3A, la courbe I 
5 (respectivement II) correspond a un indice de 
refraction n qui suit parfaitement la loi choisie 
(respectivement a un indice de refraction n trouv6 par 
la m^thode des nceuds) . 

Sur la figure 3B, la courbe I 
10 (respectivement II) correspond a un coefficient 
d' extinction k qui suit parfaitement la loi choisie 
(respectivement a un coefficient d' extinction k trouve 
par la m^thode des noeuds) . 

On vient de d^crire des exemples de 
15 1' invention. On remarquera d'une maniere plus g6n6rale 
que, dans cette derni^re, on considSre un ensemble X de 
valeurs, avec X= {ni^a, w/ ni / ..ji m/ ki # k 2/ -. / ki / .,. / k m/ 8} / ou 

ni est la valeur de l 1 indice de refraction 
(r€el) au noeud correspondant a X±, ie{l...m}, m €tant le 
20 nombre de noeuds 

ki est la valeur du coefficient d 1 absorption 
au noeud correspondant a Xi, ie {l—m}, 

e est l f £paisseur de la couche £tudi§e. 

Dans ce cas, 1 ! operation de minimisation 
25 d'erreur Er (\P, W ) revient a trouver 1' ensemble ou 
«: vecteur » X tel que Er soit minimum. 

Lorsque l'on n' impose pas de contrainte 
particulidre, la minimisation est une minimisation a 
2xm+l paramStres. On peut bien sur a j outer des 
30 contraintes afin de diminuer le nombre de variables. 
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En particulier, si l'on salt que le 
mat^riau est non absorbant, on impose ki=0 pour tout i 
de {l~m} et X devient : X={n a , n 2 .« , n A , xx^ , e} . 

Si, par tine raesure complement aire (par 
5 exemple une mesure de goniom^trie ou une mesure directe 
non optique) , on connait l'epaisseur de la couche 
consider€e avec une precision suffisante, l'^paisseur e 
n'est plus une variable et l'on a : 
X-{ni , n 2 , ... , ni , .~xi m , k x , k 2 , ... , ki , , km} . 
10 Bien entendu, les deux options pr6c6dentes 

peuvent §tre combin^es. 

On explique dans ce qui suit un mode de 
mise en oeuvre de 1' invention sous une forme plus 
g£n£rale que les exemples donnes precedemment . 
15 Soit a (X) une fonction de la longueur 

d'onde X de la lumi£re utilis§e pour les mesures. On 
peut par exemple choisir : 

a(X)=X (voir les figures 3A et 3B oil le 
spectre est uniform^ment £chantillonn6 en X) 
20 a(X)~l/X (voir la figure 2 oil le spectre 

est uniformement echantillonne en X/X) 

a{X)=ho/X ou h est la constante de Planck 
et c la vitesse de la lumi^re dans le vide, a(X) £tant 
alors homogdne k une 6nergie. 
25 Soit A l'intervalle du spectre de mesure, B 

l'intervalle spectral d€crit par les « noeuds » et C 
l'intervalle d'int6r8t. 

L'intervalle C est inclus dans l'intervalle 
B ou 6gal k cet intervalle B. De mime, l'intervalle B 
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est inclus dans l'intervalle A ou egal a cet intervalle 
A. 

On precise que chacun des intervalles A, B 
et C est du type [Om, <Xm] ou 0^ est inferieur a (Xm et il 
5 existe deux longueurs d'ondes X* et %i telles que 
OCm^at^c) cxM=a(Xi) . A titre purement indicatif et 

nullement limitatif, la figure 4 montre un exeraple des 
intervalles A, B et C et des courbes des variations de 
¥ et de ¥ en fonction de a(X) , *F £tant en fait une 

10 fonction de n*(a(A0) • 

Les cercles N repr^sentent les noeuds. 
n* est l'indice complexe, qui est exprim£ 
ici en fonction de a (A,) et dont les parties r§elle et 
imaginaire sont respect ivement not€es n(a(X)) et k 
15 (a(A,)). 

On voit ^galement, a titre purement 
indicatif et nullement limitatif, un exemple de la 
courbe des variations de n (respect ivement k) en 
fonction de a (A,) , passant par des points de coordonn^es 

20 (ni/ a±) (respect ivement (ki, a*.) , ou tx±=(X(X±) , ni=n(ai), 
ki«k(oti), l^i<m (m entier nature 1 non nul) . 

Sur la figure 4, on notera la 
correspondance entre les nceuds N, les points (ni, (Xi) et 
les points (ki, oti) . 

25 Dans le mode de mise en oeuvre consid€r£, on 

utilise un algorithme comprenant les Stapes suivantes : 

1- on fait les mesures Y sur l'intervalle A 
et 1'on choisit les m£thodes de calcul associ^es aux 
mesures (calculs ellipsom^triques ou 

30 r§f lectom^triques) ; 
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2. on choisit m nombres Oi (constituant m 
valeurs initiales de la fonction a) , i appartenant & 
{l,~.,m}, avec m>l, et {oci} c= A (leB a* correspondant aux 
« noeuds ») ; lorsque m>l, B est defini comme 1' ensemble 

5 des points a tels que min(Oi) ^0t£ max((Xi) ; lorsque 111=1, 
on a B=A ; 

3 . on choisit m valeurs initiales d' indice 
complexe n* aux m points OL± f i appartenant k {l /w ,m} ; 

4. si xc&l on choisit une loi 
10 d' interpolation qui permet de calculer 1' indice de 

refraction n(a) sur l'intervalle B a partir des points 
(Oi, ni) , i appartenant a {l,.~,m} ; si m«l, alors 
n(a)=ni(ai) sur tout l'intervalle B ; 

5. on choisit M param&tres variables avec 
15 M£2m+1 ; ces paramdtres peuvent §tre par exemple : 

les parties reelles des indices de 
refraction aux points a±,i appartenant a {l,.~,m}, ou 

- les parties imaginaires des indices de 
refraction en ces raSmes points , ou 
20 - l'6paisseur du materiau dont on cherche 

1' indice de refraction ; 

6. on choisit une fonction d'erreur 
Er^, *F) qui caracterise l'ecart entre un spectre 
mesure et un spectre theorique ; en general, l'erreur 

25 est mesuree sur l'intervalle C ; 

7. £ 1'aide d'une fonction de minimisation 
de Er(>F,\P) & M param£tres, on effectue la serie 
d'etapefe suivantes : 
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a) k 1 ' aide de la f onction d' interpolation 
des (Oti, ni) sur B, on deduit n(a) avec a appartenant k 
B ; 

b) k l'aide de n{a) de l'epaisseur 8 et des 
5 methodes de calcul de spectre, on calcule le spectre 

theorique Y (n(a),8) ; 

c) on compare et W k l'aide de 

Er=Er (*F,*F) ; si Er est suffisamment petite (c'est-t- 
dire si Er est inf 6rieure i. une valeur pr£d#f inie e) , 
10 ou bien si Er est minimale, on passe k l'€tape e) , 
sinon on passe k 1'^tape d) ; 

d) on fait varier les M paramStres de facon 
k tendre vers le minimum de Er et 1'on revient k 
l'etape a) ; 

15 e) si Er est inf€rieure k e on obtient done 

un ensemble de M paramitres tel que Er(T, T (n(a,M),e)) 
est minimum et le calcul de l'indice de refraction est 
acheve : cet indice est pris egal k celui qui a €te 
obtenu en dernier ; et si Er est sup§rieure ou egale k 

20 e, on passe l'etape 8) ; 

8. on augmente m et 1'on revient k l'etape 

2) . 

II convient de noter que la presente 
invention peut utiliser non seulement des 
25 e chant illonnages en longueur d'onde (X) mais encore des 
€chantillonnages en frequence {c/X ou, plus simplement 
lA) , en energie (he A) et, plus gin^ralement , en 
paramdtre f onction de la longueur d'onde. 

II convient en outre de noter qu'une §tape 
30 essentielle de l'algorithme (£tape 8) raentionn6e plus 
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haut) n'est pas limitee & 1'ajout d'un ensemble de 
noeuds aux noeuds d§j^ existants : elle englobe le cas 
plus general ou le nombre de noeuds augmente. Cela veut 
dire, dans un mode de realisation particulier de 
5 1' invention, qu'aprds minimisation avec 3 noeuds, si 
1'on veut passer & 6 noeuds au total, la position des 3 
anciens noeuds est « effac€e » de facon, par exemple, 3. 
avoir une densite de noeuds constante sur le spectre. En 
pratique, c'est la meilleure option, L' information sur 

10 1'ancienne position des anciens noeuds n'est pas perdue 
car la valeur des indices aux anciens noeuds est 
utilisee pour calculer la valeurs des 6 nouveaux noeuds 
(en fait 3 nouveaux noeuds plus trois anciens noeuds) . 

Ainsi, selon un mode de mise en oeuvre 

15 particulier du proceed objet de 1' invention, on peut 
augmenter le nombre de valeurs initiales de la fonction 
a en ajoutant une ou plusieurs valeurs aux valeurs 
initiales existantes ; mais, selon un mode de mise en 
oeuvre pr£f6re, on peut augmenter le nombre de valeurs 

20 initiales de la fonction a en remplacant les valeurs 
initiales existantes par de nouvelles valeurs initiales 
dont le nombre est superieur au nombre des valeurs 
initiales existantes. 

La pr£sente invention n'est pas limitee & 

25 la caracterisation de couches minces. Elle s' applique 
aussi k la caracterisation de couches epaisses. 

En outre, la prisente invention n'est pas 
limited & la caracterisation d'une seule couche, formee 
sur un substrat. Elle s' applique aussi k la 

3 0 caracterisation de deux, ou plus de deux, couches 
f ormees sur- un substrat . 
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REVEND X CATIONS 
1. Proc€d§ de caract^risation optique d'au 
moins une couche d'un materiau dans un intervalle A de 
valeurs prises par une fonction a d'une longueur d'onde 
5 optique X, lorsque X varie dans un intervalle de 
longueurs d'ondes, cette couche etant form£e sur un 
substrat, ce proce"d6 e"tant caractexise" en ce qu'il 
comprend les Stapes suivantes : 

1) on effectue un ensemble de mesures de r§f lectom^trie 
10 et/ou d'*ellipsom6trie sur 1' intervalle A, cet ensemble 

de mesures conduisant a un spectre mesure\ note" et 
1'on choisit les mSthodes de calcul associ£es a la 
nature des mesures et au type de couche a caract€riser, 

2) on choisit m valeurs initiales <Xi .-a* de la fonction 
15 a, appartenant a cet intervalle A, m Stant un nombre 

entier au moins egal a 1, et l'on d£finit un 
intervalle B comme £tant 1' ensemble des points a de 
1' intervalle allant du plus petit au plus grand des 
nombre s CCi-O^, lorsque m est superieur a 1, et comme 
20 etant 1' intervalle A lorsque m vaut 1, 

3) on choisit m valeurs initiales complexes d'un indice 
de refraction complexe n*=n+jk aux m points Ok, i allant 
de 1 a m, 

4) lorsque m est different de 1, on choisit une loi 
25 d' interpolation qui permet de calculer 1> indice de 

refraction n(a) du materiau sur 1' intervalle B, a 
partir des points (a±, n ± ) , avec ni^nfc^), i allant de 1 
a m, et lorsque m vaut 1, n(a) est pris egal au nombre 
ni(Oi) sur tout 1' intervalle B, 
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5) on choisit M parametres variables, M £tant infer ieur 
ou egal a 2m+l, 

6) on choisit une fonction d'erreur Er(\F, W) qui 
caracterise l'ecart entre un spectre mesure et un 

5 spectre theorique *P, 

7) k l'aide d'une fonction de minimisation de Er (*F, W) 
k M parametres, on effectue la serie d'€tapes 
suivante : 

a) k l'aide de la loi d' interpolation des (a^ni) sur 
10 l'intervalle B, on deduit n(a) , a appartenant k B, 

b) k l'aide de n(a) et de l'epaisseur 8 de la couche, 
et de m^thodes de calcul de spectres, on calcule un 
spectre theorique *F (n(a) ,e) , 

c) on compare ¥ et ¥ k l'aide de ErPF, W) et, si 
15 Er(¥, *P) est suffisamment petite, c'est-a~dire 

inferieure k une valeur pr^definie e, ou est minimale, 
on passe k l'6tape e) , sinon on passe k 1'etape d) , 

d) on fait varier les M parametres variables de f aeon k 
tendre vers le minimum de Er C¥, *F) et 1'on retoume k 

20 1' §tape a) , 

e) si Er(¥, W) est inferieure k e, on obtient done un 
ensemble de M parametres variables, pour lequel Er(\P, 
¥(n(oc,M),8)) est minimum et l'indice de refraction est 
alors pris egal k celui qui a 6t€ obtenu en dernier, et 

25 si Er(¥, *F) est superieure ou egale k e on passe k 
l'€tape 8) , 

8) on augmente le nombre m de valeurs initiales de la 
fbnption a et l'on retoume a l'^tape 2). 
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2. Proc6d€ selon la revendication 1, dans 
lequel on augmente le norabre de valeurs initiales de la 
fonction a en ajoutant une ou plusieurs valeurs aux 
valeurs initiales existantes. 
5 3. Proc6d€ selon la revendication 1, dans 

lequel on augmente le nombre de valeurs initiales de la 
fonction a en remplagant les valeurs initiales 
existantes par de nouvelles valeurs initiales dont le 
nombre est supirieur au norabre des valeurs initiales 
10 existantes. 

4. Proc£d6 selon l'une quelconque des 
revendications 1 a 3, dans lequel chaque loi 
d 1 interpolation est choisie parmi les lois 
d 1 interpolation lin^aires , les lois d 1 interpolation 

15 cubiques, les lois d» interpolation polyndmiales et les 
lois d 1 interpolation par exemple de type fonction 
spline. 

5. Procede selon 1'une quelconque des 
revendications 1 k 4 , dans lequel les valeurs 

20 initiales de la fonction a sont reparties de mani^re 
uniforme dans 1' intervalle A, la repartition des noeuds 
etant ainsi homog^ne. 

6. Proc€d§ selon l'une quelconque des 
revendications 1 & 5, dans lequel a(X) est choisi parmi 

25 X, 1/X et hcA, ou h est la constante de Planck et c 
est la Vitesse de la lumiSre dans le vide. 

7. Proc£d6 selon l'une quelconque des 
revendications 1 k 6, dans lequel on mesure l'erreur, It 
l'gtape 6),' sur un intervalle d'int^ret C qui est 

30 inclus dans 1'intervalle B ou ggal £ cet intervalle B. 
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8. Precede selon l'une quelconque des 
revendications 1 & 7, dans lequel les M param^tres 
variables sont les parties replies des indices de 
refraction aux points (Xi, i allant de 1 am. 
5 9. Precede selon l'une quelconque des 

revendications 1 £ 7, dans lequel les M paramStres 
variables sont les parties imaginaires des indices de 
refraction aux points a*,, i allant de 1 £. m. 

10. Procede selon l'une quelconque des 
10 revendications 1 k 7, dans lequel les M paramStres 

variables sont constitue*s par 1'epaisseur du materiau 
dont on cherche 1' indice de refraction. 

11. Procede de caracterisation optique d'au 
moins une couche d'un materiau dans un intervalle de 

15 longueurs d'ondes [X min, X max], cette couche etant 
formee sur un substrat, ce procede etant caracterise en 
ce que : 

- on effectue un ensemble de mesures de r£f lectometrie 
et/ou d'ellipsometrie, cet ensemble de mesures 

20 conduisant Tin spectre mesure, note *F, 

- on choisit m longueurs d'ondes initiales X x Jk^ 
appartenant k cet intervalle, m etant un nombre entier 
au moins egal cl 1, on associe, si chaque longueur 
d'onde, un indice de refraction, 

25 - on choisit une loi d 1 interpolation au moins pour 
I'indice de refraction du materiau, pour les longueurs 
d'ondes comprises entre les longueurs d'ondes 
initiales Xi-Xn, 

- on choisit M paramStres initiaux, M etant au moins 
30 egal & m, k savoir un indice de refraction initial n A 
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pour chaque longueur d'onde initiale X±, 1 ss i <; m, 
les longueurs d'ondes initiales 6tant choiaies de 
mani^re a. pouvoir determiner par interpolation au 
moine 1 1 indice de refraction pour toute longueur 
5 d'onde de I'intervalle [X min, X max], les couples 
(Xi/ni) etant appeles nceuds, 

— on choisit des mithodes de calcul de re"f lectom£trie 
et d 1 ellipsom£trie, 

— on choisit aussi une fonction d'erreur Er, 
10 representative de l'ecart entre deux spectres et 

*F 2 , les spectres *Fi et *F 2 §tant calculus ou raesure*s 
sur un nombre de points sup£rieur au nombre m de 
noeuds, 

— a I 1 aide des m longueurs d»ondes initiales, des M 
15 param&tres initiaux et de la loi d ! interpolation, on 

met en oeuvre le processus d 1 optimisation suivant : 

— on determine un spectre th§orique, note 

T , en fonction des mithodes de calcul choisies, et de 
1' indice d£duit par interpolation de sa valeur en A*, i 
20 allant de 1 a m, sur le spectre [tanin, Amax] , 

— on determine l'erreur Er (*F, T) , entre le spectre 
mesur€ et le spectre theorique 

— on minimise cette erreur en faisant varier la 
position des valeurs des indices inconnus et/ou 

25 l'epaisseur de la couche et/ou les valeurs des indices 
de refraction aux longueurs d'ondes initiales, et l'on 
obtient un spectre, 
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— on ajoute des longueurs d 1 ondes aux longueurs d 1 ondes 
initiales A* ~ 7^, lea longueurs d' ondes ajout€es 
const ituant de nouveaux noeuds, 

— on r^pdte le proc£d6 en choisissant un nombre m' de 
5 longueurs d' ondes initiales, m' €tant sup€rieur a m, 

et M' paramdtres initiaux, M' £tant sup#rieur a M, 
jusqu'a ce que la precision sur chaque spectre ainsi 
repr6sent£ au mieux soit 6gale a une precision 
pr€d€f inie. 

10 12. Proc£d£ selon la revendication 11 , dans 

lequel m est au moins §gal £ 2. 

13. Proced£ selon la revendication 11, dans 
lequel m est €gal a 1 et l'on choisit des indices de 
refraction initiaux €gaux. 

15 14. Proc£d6 selon l'une quelconque des 

revendications 11 a 13, dans lequel le matSriau est non 
absorbant et le notnbre M est £gal a m, le coefficient 
d' extinction du materiau etant pris €gal a 0. 

15. Proc§de selon I 1 une quelconque des 

20 revendications 11 a 13, dans lequel : 

— M est au moins €gal a 2 m, 

— on choisit en outre une loi d 1 interpolation pour le 
coefficient d 1 extinction du materiau, 

— pour chaque longueur d'onde initiale A*, 1 * i s m, 
25 on choisit en outre un coefficient d' extinction initial 

k if les longueurs d' ondes initiales 6tant en outre 
choisies de manidre a pouvoir determiner par 
interpolation le coefficient d 1 extinction pour toute 
longueur d'onde de 1 1 interval le [A, min, X max], 
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- dans le processus d 1 optimisation, on minimise 
l'erreur en faisant varier en outre les valeurs des 
coefficients d ! extinction aux longueurs d'ondes 
initiales et les longueurs d'ondes ajoutees sont en 
5 outre plac^es de fagon a representer au mieux le 
spectre du coefficient d 1 extinction du materiau. 

16. Procede selon la revendication 15, dans 
lequel m est egal a 1 et l'on choisit des indices de 
refraction initiaux egaux et des coefficients 

10 d' extinction initiaux egaux. 

17. Procede selon l'une quelconque des 
revendications 11 a 16, dans lequel la couche de 
materiau est mince, c'est-a-dire a une epaisseur 
inferieure a la longueur de coherence de la lumi^re 

15 utilis^e pour les mesures, on choisit un param^tre 
initial supplement air e , a savoir une Epaisseur initiale 
de couche, et dans le processus d 1 optimisation on 
minimise l'erreur en faisant en outre varier la valeur 
de l 1 Epaisseur de couche. 

20 18. Procede selon I 1 une quelconque des 

revendications 11 a 16, dans lequel la couche de 
materiau est epaisse, c'est-a-dire n'est pas mince, et 
M est au plus egal a 2 m. 

19. Precede selon l'une quelconque des 
25 revendications 11 a 16, dans lequel 1" epaisseur de la 

couche de materiau est connue avec une precision 
predefinie et M est au plus egal a 2 ra. 

20. Procede selon l'une quelconque des 
revendications 11 a 19, dans lequel chaque loi 

30 d 1 interpolation est choisie parmi les lois 
d 1 interpolation lineaires, les lois d» interpolation 
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cubiques, les lois d 1 interpolation polyn6Ttiiales et les 
lois d ! interpolation de type fonction spline. 

21. Proc£d6 selon l'une quelconque des 
revendi cat ions 11 & 20 , dans lequel la repartition des 
5 noeuds est homogene. 
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